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We revisit the problem of the dispersion of a single hole injected into a quantum antiferromagnet. 
We applied a spin-density-wave formalism extended to large number of orbitals and obtained an 
integral equation for the full quasiparticle Green's function in the self-consistent "non-crossing" 
Born approximation. We found that for t/J ^ 1, the bare fermionic dispersion is completely 
overshadowed by the self-energy corrections. In this case, the quasiparticle Green's function contains 
a broad incoherent continuum which extends over a frequency range of ~ 6t. In addition, there exists 
a narrow region of width 0{JS) below the top of the valence band, where the excitations are mostly 
coherent, though with a small quasiparticle residue Z ~ J/t. The top of the valence band is located 
at {it/2, tt/2). We found that the form of the fermionic dispersion, and, in particular, the ratio of the 
effective masses near {it/2, tt/2) strongly depend on the assumptions one makes for the form of the 
magnon propagator. We argue in this paper that two-magnon Raman scattering as well as neutron 
scattering experiments strongly suggest that the zone boundary magnons are not free particles since 
a substantial portion of their spectral weight is transferred into an incoherent background. We 
modeled this effect by introducing a cutoff in the integration over magnon momenta. We found 
analytically that for small Qc, the strong coupling solution for the Green's function is universal, and 
both effective masses are equal to {4JS)^^. We further computed the full fermionic dispersion for 
J/t = 0.4 relevant for Sr2Cu02Cl2, and t' — — 0.4J and found not only that the masses are both 
equal to (2J)~^, but also that the energies at (0,0) and (0, tt) are equal, the energy at (0,7r/2) is 
about half of that at (0, 0), and the bandwidth for the coherent excitations is around 3 J. All of these 
results are in full agreement with the experimental data. Finally, we found that weakly damped 
excitations only exist in a narrow range around (7r/2,7r/2). Away from the vicinity of {n /2,tt/2), 
the excitations are overdamped, and the spectral function possesses a broad maximum rather than 
a sharp quasiparticle peak. This last feature was also reported in photoemission experiments. 



I. INTRODUCTION 

The dispersion of a single hole in a quantum an- 
tiferromagnet is one of the issues in the field of 
high-temperature superconductivity which has attracted 
a substantial amount of interest over a number of 
years jl|~[ig|]. The parent compounds of the high-Tc 
materials are quantum Heisenberg antiferromagnets as 
was demonstrated by neutron scattering NMR ||2|] 
and Raman experiments. The antifcrromagnctic 
spin ordering strongly modifies the electronic dispersion 
which by all accounts is very different from what one 
would expect from band theory calculations. Upon hole 
doping, short-range antiferromagnetism gradually disap- 
pears, and the overdoped cuprates possess an electronic 
dispersion which is consistent with band theory predic- 
tions . How the electronic spectrum evolves with dop- 
ing is currently a subject of intensive experimental and 
theoretical studies [p6|-p9[. As an important input for 
these studies, one needs to know what happens in the 
limit of zero doping when a single hole is injected into a 
quantum antiferromagnet. 

The dispersion of a single hole in an antiferromagnet 
has been intensively studied experimentally and theoreti- 
cally . Experimental information comes from photoemis- 
sion experiments on the half-filled Sr2Cu02Cl2 which 
is not a high-Tc superconductor, but contains the same 



Cu02 planes as the high-Tc materials 31|. Most of 
the theoretical analysis was performed in the framework 
of the t — J and Hubbard models which are widely be- 
lieved to adequately describe the low-energy physics of 
the underlying three-band model p]" p^ , p^ -pj|] . Early 
analytical and numerical computations were performed 
in the antiferromagnetically ordered phase and for the 
case when a hopping is only possible between nearest 
neighbors |^,|^,|l^,|l3| . These studies have shown that 
in the strong coupling limit (large U limit in the Hubbard 
model or t ^ J limit in the t — J model), the Green's 
function of a single hole has the form 



G{k,u) 



Z 



Eh 



Ginc{k,Uj) , 



(1) 



where the coherent part is confined to scales smaller than 
2 J, while the incoherent background stretches upto a few 
t. The quasi-particle residue of the coherent piece is small 
and scales as Z oc J/t in the limit t ^ J. The disper- 
sion Ek has a maximum at fc = (7r/2,7r/2) and symme- 
try related points. All calculations have demonstrated 
that the dispersion around this point is very anisotropic 
with a substantially larger mass along the (0, tt) to (tt, 0) 
direction than along the Brillouin zone diagonal. For 
t/J~ 2.5 relevant to cuprates, the ratio of the masses 
is about 5 — 7 in the t — J model (without a three-site 
term) |^ , and it is even larger in the Hubbard model due 
to the presence of the bare dispersion J{coskx + cosky)^ 



1 



which yields an extra contribution to the mass along the 
zone diagonal ||2C[] . 

It turns out, however, that the experimental results for 
Sr2Cu02Cl2 10,0 arc rather different from these pre- 
dictions. Although the photoemission data have demon- 
strated that the maximum of Ek is at A; = (7r/2,7r/2) 
consistent with the theory, the experimentally measured 
ratio of the masses is close to one in clear disagreement 
with the theoretical predictions. Moreover, the data show 
that the coherent peak in the spectral function exists only 
in a narrow region around (7r/2,7r/2) while away from 
this region, the hole spectral function is nearly feature- 
less. This implies that the fermionic excitations become 
overdamped already at energies which are substantially 
smaller than 2 J. 

After the data were reported, several attempts have 
been made to improve the agreement between theory 
and experiment. One scenario was put forward by re- 
searchers working on the "gauge theory" approach to 
cuprates most recently by Laughhn [|3j. He ar- 

gued that the isotropy of the dispersion together with the 
observed mostly incoherent nature of the electronic exci- 
tations are signatures of a spin-charge separation. For a 
state where spin and charge degrees of freedom are de- 
scribed by separate quasiparticles (spinous and holons, 
respectively), the electron Green's function is just a con- 
volution of the spinon and holon propagators. It does 
not have a pole which normally would be associated with 
the coherent part of G{k,uj), but rather a branch cut 
which describes fully incoherent excitations. Laughlin 
argued that since spinon and holon energies are well sep- 
arated (the spinon energy has an overall scale of J, while 
the holon energy is 0{t)), the position of the branch cut 
virtually coincides with the spinon dispersion. In the 
mean-field theory for the spin-charge separated state, the 
spinon energy has the form 

Os2fc,)l/2 , 
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where Csw 1.6J is the spin-wave velocity in a 2D S — 
1/2 antiferromagnet. This dispersion has an isotropic 
maximum at k = (7r/2,7r/2), a bandwidth of 2. 2 J and 
equal energies for k = (0, 0) and (0, tt) - all of these fea- 
tures are consistent with the data together with the near 
absence of the quasiparticle peak. 

An obvious weakness of the mean-field analysis of 
spinous and holons is that it neglects the effects due to 
a gauge field. Beyond the mean-field level, a gauge field 
may glue spinous and holons into a bound state thus ren- 
dering the electron as a coherent quasiparticle. Laughlin 
conjectured that the confinement takes place only below 
T/v, while the experimental data were actually collected 
at T = 350K which is lOOK above the Neel temperature. 
He then proposed that if measurements are done at much 
lower temperatures, they should yield an anisotropic dis- 
persion consistent with the results obtained in the or- 
dered state with no spin-charge separation. 

Another, more conventional approach to the single hole 
problem assumes that there is no spin-charge separation 



at any T, and that the experimental data in fact reflect 
the behavior of the hole dispersion in the antiferromag- 
nctically ordered phase ]l5|-p^. Within this approach, 
the discrepancy with the data is mainly attributed to the 
fact that the original model did not contain a hopping 
term t' between next-nearest neighbors (and, possibly, 
also between further neighbors). The presence of the a 
finite t' term in the Hubbard model is justified, at least 
partly, by studies which derived an effective one-band 
model from the underlying three-band model by com- 
paring the energy levels around the charge transfer gap 
iQ. These studies predicted that the second- neighbor 
hopping is about t' — —0.2t. By itself, this hopping is 
small compared to t. However, in an antiferromagnetic 
background, the hole can only move within the same sub- 
lattice, otherwise the antiferromagnetic ordering is dis- 
turbed. The hopping term t' connects the sites within 
the same sublattice, and therefore is not affected by anti- 
ferromagnetism. On the contrary, the t term contributes 
to the hopping within a sublattice only via the creation 
of a virtual doubly occupied state which costs the energy 
U. As a result, the part of the dispersion is rescaled 
and becomes of order t"^ /U = 0{J). One therefore has to 
compare t' not with t but rather with J. For J/t ^ 0.4, 
we then obtain t' = —0.5 J, which immediately implies 
that the corrections due to t' are actually quite relevant. 

It has been mentioned several times in the literature 
that the inclusion of i' = —0.5 J into the Hubbard model 
yields a good agreement with the experimental data al- 
ready at the mean- field level |l|,|l6|. Indeed, the mean- 
field spin-density-wave (SDW) formula for the hole dis- 
persion at large U is 

Ek = — J(cos kx + cos kyf' — 4i' cos k^ cos ky . (3) 

For t' = —0.5 J, this formula transforms into 

Ek = - J(cos^ kx + cos^ ky) (4) 

(here we assumed that the chemical potential is at the top 
of the valence band). This dispersion possesses two equal 
effective masses if one expands around the maximum at 
(7r/2,7r/2), and has a a local maximum at (0,7r/2) with 
E = — J. Both of these results are consistent with the 
most recent data by LaRosa et al. BL] . Furthermore, the 
data show that the energies at (0,0) and (0,7r) are both 
equal to —2J. This also agrees with the photoemission 
data (30|j3^. 

The conventional mean-field SDW-type approach also 
possesses the weakness that it predicts fully coherent ex- 
citations upto 2 J. The data, however, demonstrate that 
away from the vicinity of (7r/2, 7r/2), the coherent part of 
the dispersion is almost completely overshadowed by the 
incoherent background. Earlier studies [ p^ which went 
beyond the mean-field level have demonstrated that self- 
energy corrections reduce the quasiparticle residue thus 
transferring part of the spectral weight into the inco- 
herent background. However, these corrections also ef- 
fectively decrease t' and thus render the spectrum more 
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anisotropic (see Fig. ^ and below). From this per- 
spective, the observed isotropy of the dispersion around 
(7r/2, 7r/2) is attributed in a conventional approach to 
some fine tuning of both J /t and t' / J and is therefore 
completely accidental 

In this paper we show that in a certain limit speci- 
fied below, the near-degeneracy of the spectrum around 
(7r/2, 7r/2) turns out to be a fundamental, universal prop- 
erty of a single hole in an antiferromagnet, indepen- 
dent of the details of the physics at atomic scales. Our 
key point is this: in all previous studies which yielded 
anisotropic spectra, it was assumed that magnons behave 
as free particles for all momenta. In this case, the inte- 
gral over the magnon momenta in the self-energy term 
runs over the whole magnetic Brillouin zone (MBZ). On 
the other hand, Raman studies of the two-magnon profile 
in the insulating parent compounds of high Tc materials 
have demonstrated that the width of the two-magnon 
peak is much broader than one would expect for free 
magnons 1 23 j2^ . The dominant contribution to this peak 
comes from the magnons near the boundary of the MBZ. 
Complimenting these findings, neutron scattering exper- 
iments on La2CuOi ||3^ have shown that about half of 
the spectral weight of the quasiparticle peak for the zone 
boundary magnons is transferred into a broad incoherent 
background. 

It has been suggested that the broadening is due 
to the strong interaction between these magnons and 
phonons p7| , p8[ . This interaction is finite and not nec- 
essary small at T = contrary to the magnon-magnon 
interaction which gives rise to an incoherent part of the 
magnon spectral function only at finite T and is irrele- 
vant for T< J 19). 

In this situation, it seems reasonable to assume that 
the contribution from the zone boundary magnons to the 
electronic self-energy is substantially reduced compared 
to what one would obtain for free spin waves. This how- 
ever is true only for zone-boundary magnons. For long- 
wavelength magnons, the magnon-phonon vertex scales 
linearly with the magnon momentum, and the incoher- 
ent part of the magnon propagator is small. The simplest 
way to model this effect is to introduce an upper cutoff 

in the integration over magnon momenta. Naively, 
one might expect that the hole dispersion would strongly 
depend on qc- However, we will demonstrate that at 
large t/ J, when the bare dispersion is irrelevant, only the 
quasiparticle residue does depend on while the effec- 
tive masses are in fact independent of qc in the limit when 
qc is sufficiently small. We explicitly show that in this 
limit, both masses turn out to be equal to 1/2 J. The dis- 
persion near (7r/2,7r/2) is then isotropic and has a form 
Ek = —Jk^ where k is the deviation from (7r/2,7r/2). 
Furthermore, we show that for a certain range of qc the 
inclusion of t' = —0.5 J extends the region where the two 
masses are approximately equal to basically all values of 
t/J. This last result allows us to correctly reproduce the 
measured hole dispersion in Sr2Cu02Cl2- 



The paper is organized as follows. In the next sec- 
tion, we outline the formalism and derive the integral 
equation for the quasiparticle Green's function by ex 
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panding around the mean-field SDW solution. In Sec 
we present our analytical results in the large t/J limit. 
In this section we also discuss the role of the vertex 
corrections to the spin-fermion vertex. In Sec. we 
present the results of the numerical solution of the self- 
consistency equation for the quasiparticle Green's func- 
tion for different values of J/t. Sec. ^ contains a sum- 
mary of our results. 



II. THE FORMALISM 

As mentioned in the introduction, our starting point 
for the description of the insulating parent compounds of 
the high- Tc materials is the effective 2D one-band Hub- 
bard model H-pl, given by 



(5) 



Here a is the spin index and ti^j is the hopping integral 
which we assume to act between nearest and next-nearest 
neighbors (t and t' , respectively). Throughout the paper 
we assume that the ground state of the Hubbard model 
is antiferromagnetically ordered. In this situation, a way 
to calculate the spectral function in a systematic per- 
turbative expansion is to extend the Hubbard model to 
a large number of orbitals, n = 2S, and use a l/S* ex- 
pansion around the mean-field SDW state |4^. The 1/5* 
expansion for the Hubbard model has been discussed sev- 
eral times in the literature and we will use it here 
without further clarification. To obtain the mean-field 
solution, one introduces an antiferromagnetic long range 
order parameter Sz = (cj^.C/t+g) and uses it to decou- 
ple the interaction term in Eq.(^. Diagonalizing then 
the quadratic form by means of a unitary transformation 
one obtains two electronic bands for the conduction and 
valence fermions, whose energy dispersion is given by 



where 



± _ £*; ± ffc+C 



-4i5(cos kx + cos ky) 



A2 



A = U{Sz 



8t'S cos kx cos ky — fi 



(6) 



(7) 



Here -E^'^ is the dispersion of the conduction and va- 
lence fermions, respectively, ek is the dispersion of free 
fermions, fj, is the chemical potential, and (Sz) is the 
sublattice magnetization. To facilitate the 1/S expan- 
sion, we also introduced t = t/2S and P = t' /2S. In the 
large-J7 limit which we only consider, one can expand the 
square root and obtains 
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FIG. 1. The lowest order self-energy correction for the va- 
lence fermions in the SDW model. The solid and dashed lines 
are the bare propagators of conduction and valence fermions, 
respectively. The wavy line describes transverse spin fluctua- 
tions. 



f, — ±A±2JS{coskx + cosky) 



-8t'S cos kr cos L, 



(8) 



where J = 4P/U. At half-filling, the chemical potential 
can be set to the top of the valence band (/i = — A); for 
S = 1/2 we then reproduce Eq.(^ from the introduction. 

At infinite S, the mean-field approach is exact. At 
finite S, the bare Green's function is rcnormalized due 
to the interaction with spin waves. The lowest order 
self-energy corrections for valence fermions are given by 
the diagrams in Fig. |^. The solid and dashed lines in 
these diagrams are the propagators of conduction and 
valence fermions, respectively. The wavy lines describe 
transverse spin fluctuations which in the SDW approach 
are collective modes of electrons. These collective modes 
correspond to the poles of the transverse susceptibility, 
and are obtained by summing up an infinite RPA se- 
ries in the particle-hole channel with the total momen- 
tum equal to either zero or Q. The interaction vertices 
between fermionic quasiparticle and magnons have been 
calculated previously |Q. In the strong coupling limit 
they are given by 



^cc,vvik,q)= 



k+q )''h 



X-) 



'i + q)Vq 



where rjq and fjq are given by 
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Vq 
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, (10) 



and Vq = {cosqx + cosgy)/2. 

We see that there are two types of vertices: ^cv.vc 
which describes the interaction between conduction and 
valence fermions, and ^cc,vv which involves either only 
valence or only conduction fermions. Apparently, the 
second diagram in Fig. ^ is more relevant since the ver- 
tex which involves both conduction and valence fermions 
scales as U. However, incident and intermediate fermions 
in this diagram belong to different bands and are there- 
fore separated by a large, momentum independent gap 
A ~ US. As a result, the first diagram mostly con- 
tributes to the gap renormalization, which is exactly can- 
celled by a renormalization of (Sz) such that the fully 
rcnormalized gap equals 2U S as it indeed should be for 
the large U Hubbard model Expanding this di- 

agram in J/C/, we also obtain a momentum dependent 




FIG. 2. The lowest order vertex correction for the vertex 
between fermions and transverse spin fluctuations. The dia- 
gram with only one wavy line is absent in the ordered state 
as it does not conserve the spin. 



term of 0{J) which contributes a regular 1/S correction 
to the bare dispersion. 

The first diagram in Fig. |^ involves only valence 
fermions. Here the vertex is reduced from U due to the 
coherence factors and scales as t. At the same time, both 
incident and internal quasiparticles are only 0( J) away 
from the Fermi surface which implies that the denomina- 
tor scales as J. The total contribution from the second 
diagram then behaves as JS{t/ J^/S)^ and in addition 
is strongly momentum dependent. Since the bare dis- 
persion is of order JS, the relative self-energy correction 
from the second diagram scales as {t/J^/S)"^ and is small 
only for extremely large S. For physically relevant values 
of the spin, the expansion parameter is obviously large, 
and one certainly cannot restrict with the second order 
in perturbation theory. 

We now formulate precisely under which conditions we 
carry out the calculations. We assume that 3> 1 and 
neglect all regular corrections in 1/S. At the same time, 
we assume that t/ J^/S ^ 1 and sum up an infinite series 
of diagrams in this parameter. The restriction to large 
S allows us not only to neglect the self-energy diagrams 
which involve both valence and conduction fermions, but 
also to neglect the quantum corrections to the spin-wave 
propagator. At half-filling, these regular 1/5 corrections 
can, with good accuracy, be absorbed into the renormal- 
ization of the hopping term and the exchange interaction 
which are both input parameters for our calculations. 

The next step is to select the series of diagrams which 
have to be summed up. To lowest order in perturba- 
tion theory, both self-energy and vertex corrections are 
equally relevant: the self-energy correction yields a rel- 
ative contribution of {t/J^/S)"^, while the leading order 
vertex correction shown in Fig. |2| yields a relative factor 
{t/Jy/Sf which is even larger. This result, however, 
changes if we estimate the strength of the self-energy 
and vertex corrections in a self-consistent manner, i.e., 
by considering all internal Green's functions and all ver- 
tices in the diagrams in Figs. |l] and ^ as full ones. This 
in turn yields self-consistent equations for the full self- 
energy and the full vertex. Our self-consistent calcula- 
tion of the self-energy correction is similar to the one 
performed by Kane, Lee and Read (KLR) H. FoUow- 
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a) 



b) 



FIG. 3. a) The self-energy is given by an infinite sum of 
"non-crossing" diagrams, b) The Dyson equation which to- 
gether with the self-energy in a) yields Eq.(|ll|). 

ing KLR, we assume that the dominant pole approxi- 
mation for the full fermionic Green's function is valid 
upto energies of the order of the typical spin wave en- 
ergy, i.e., the full Green's function can be approximated 
as Z/{uj — Ek) where Ek = 0{JS) (we later confirm this 
result by explicit calculations). Substituting this form 
into the self-energy term and performing standard ma- 
nipulations we obtain for t/ J\/S ^ 1 the self-consistency 
condition PZ"^^/ J'^S ~ 1, where $ stands for the ver- 
tex renormalization. It is essential that there is only one 
power of $ in this relation as only one of the two vertices 
in the self-energy diagram gets renormalized. On the 
other hand, in the vertex correction diagram, all vertices 
should be considered as full ones, and the self-consistency 
condition yields {PZ'^^/J^S)'^^^ ^ 1. Comparing these 
two conditions, we obtain Z - Jy/S/t and $ = 0(1). 
The result for Z is consistent with the one obtained by 
KLR. Clearly then, the self-energy corrections are more 
relevant than the vertex corrections since the former re- 
duce the quasiparticle residue to a parametrically small 
value, while the latter only change the vertex by a fac- 
tor of order 0(1). Though the vertex corrections do not 
contain a factor 1/S, it seems reasonable to assume that 
they just change the overall amplitude of the vertex but 
do not introduce any new physics. We therefore first ne- 
glect all vertex corrections and obtain the full self-energy 
and thus the full Green's function in the self-consistent 
Born approximation We then use the solution for 

the full Green's function to estimate the relative strength 
of the vertex corrections. We find that the vertex correc- 
tions change the vertex by roughly 20% and therefore can 
be neglected with reasonable accuracy. 

In the Born approximation, the full self-energy is 
diagrammatically given by an infinite series of "non- 
crossing" diagrams (see Fig. |^). Summing up this series, 
we obtain that the full self-energy has the same form as in 
second-order perturbation theory, but the Green's func- 
tion for the intermediate fermion is now replaced by the 
full one. The full Green's function is then obtained from 
the Dyson equation (see Fig. ||b) and is analytically given 
by 



G-'{k,u;) ^cj- {El 

^ dn ^{k, q)G{k + g, CO + n)F{q, n) , 

where F{q, u) is the spin- wave propagator, and 



(11) 



^{k,q)^^t^32Si 



t2 



The integration over the magnon momentum runs over 
the whole MBZ. 

Eq.(pl]) is similar to the one derived earlier for the 
t ~ J model [|,|,|| with the only difference that Eq.|ll| 
contains the bare dispersion E'^ . This dispersion is indeed 
also present when one derives the t — J model from the 
Hubbard model at large U . However, it is due to the 
three-site term which is usually omitted in the effective 
t — J Hamiltonian . 

As we discussed in the introduction, the quasiparticle 
spectral weight of the short-wavelength magnons in the 
parent compounds of the high- Tc materials is likely to be 
strongly reduced as demonstrated by Raman and neutron 
scattering experiments. To account for this effect, we 
adopt a semi-phenomenological approach and introduce 
a cutoff, qc, in the integration over magnon momenta in 
the r.h.s. of Eq.(pl]). We assume that for q > q^ the 
magnon spectral weight disappears into a broad back- 
ground, and neglect the contribution to the self-energy 
from these q. On the other hand, for q < q^ we assume 
that the magnons are just free particles. Furthermore, 
for our analytical considerations, we will assume that qc 
is rather small such that we can expand the dispersion 
of fermions and the spin-fermion vertex to linear order 
in the magnon momentum. This last assumption is not 
well justified as the magnitude of qc is unknown. No- 
tice, however, that expanding upto leading order in q^ 
we obtain two equal effective masses which are universal 
and independent of qc- The smallness of qc is then only 
needed for the corrections to these universal results to be 
small. 

For free spin waves, we have 

Fi<l^^)^7^ , (12) 



n-, 



where lo, 



q — 4JS'^/1 — v"^ is the spin- wave spectrum. The 



magnon propagator has a pole in the lower half-plane of 
ri. In this half-plane, the mean-field fermionic Green's 
function G(A:, w) = {lu— (E^ — ^) + i6sgnjjj)~'^ is free from 
nonanalyticities since E^ — fi <Q. We assume, following 
KLR, that the full G{k,uj) is also analytic in the lower 
half- plane of fl. Then one can straightforwardly perform 
the integration over magnon frequency in Eq.(nl|) and 
obtain 



G-\k,Lu)=u;-iE]l-^i) 

^{k,q)G{k + q,UJ 



(13) 
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We first present our analytical results for the full Green's 
function in certain limiting cases, and then present the 
full numerical solution of Eq.(p^. 



III. ANALYTICAL RESULTS 

We obtain the analytical sol ution of Eq.(|13|) in two 
different ranges of w. In Sec 
self-consistency equation in the limit 
where LOmax is the highest frequency at which the full 
Green's function first acquires a finite imaginary part, 
and A = JS{i/ J\fSY/^ . We show that for \u!-ojmax\ > 
A the excitations are purely i ncohe rent and extend over 
a region of ~ &t^/2S. In Sec. Ill B we then consider the 
case \uj — 0Jmax\ ^ A. In this frequency range we find 
coherent excitations which exist up to energies of 0{JS) 
down from the maximal frequency. 



III A we first solve the 

^ynax\ > A, 



A. Incoherent part of the excitation spectrum 

We first observe that the interaction vertex in Eq. jl^ ) 
has an overall scale of (t^/S)'^, while the quasiparticle 
Green's function behaves as 1/uj at very large frequen- 
cies (here, and in the following, we shifted the frequency 
by the mean-field chemical potential, fi = —A). Obvi- 
ously than, for uj ^ t^/S, the perturbative expansion 
in the spin-fermion interaction is convergent, and the 
density of states (DOS) is exactly equal to zero, as in 
the mean-field theory. When uj is reduced to the scale 
of tV^, the lowest-order self-energy term ~ PS/uj be- 
comes of the same magnitude as the frequency in the 
bare Green's function, i.e, the expansion parameter is 
0(1). We show that in this frequency range there exists 
a critical value of uj below which perturbation theory be- 
comes non-convergent and there appears a finite DOS. It 
is essential that for small J\/S/i, the critical frequency 
is still much larger than the magnon frequency such that 
one can neglect ujq and compared to uj in the r.h.s. of 
Eq. (|l3|) . The self-consistency equation then reduces to a 
conventional integral equation 



4^ 



^ik,q-k)G{q,u;) (14) 



in which the dependence on the external momentum is 
only present in the interaction vertex. Furthermore, we 
assume that uj ~ t^/S is larger than the total magnon 
bandwidth, including the incoherent part. In this situa- 
tion, the integration over q runs over the whole MBZ. 

Before we present the solution of Eq. ( p^ ) , it is instruc- 
tive to consider a simplified version of this equation in 
which ^!{k,q — k), which is a smooth function of the 
fermionic momentum, is just substituted by some con- 
stant PS. The equation for the full G{uj) then reduces 
to 



G-\uj) =uj-t SG{uj) 



(15) 



where t/t — 0(1). Solving this algebraic equation, we 
obtain for positive uj 



G{uj) 



UJt, 



where uj„ 



2WS. We see that for uj > 



(16) 



ix^ fhe 

Green's function is real. This is the frequency range 
where perturbation theory is valid. For uj < ujmax, how- 
ever, the expression under the square root is negative, 
and the solution possesses a finite imaginary part which 
gives rise to a finite DOS. The total width of the DOS is 

obviously W — 2ujmax = 4t\/S'. 

We now solve Eq.(p^ with the actual ^{k, q — k). We 
introduce a new function fk{uj) via 



G^'iuj) ^ UJ Muj') 



(17) 



Substituting this into Eq.(14), we obtain 
d^q 



1 



1 - 1^: 



q~k 



(18) 



where we defined a ~ 32PS/uj'^. 

The general solution of Eq. (^^ can be obtained by ex- 
panding in the eigenfunctions of the D^h symmetry group 
of the square lattice. The solution is in general rather 
cumbersome because the vertex contains a fc— dependent 
term in the denominator. However, it is easy to verify 
that the expression in the square brackets vanishes when 
i^q-k ^ 1- The dominant contribution to the r.h.s. of 
Eq.(p^ then comes from the region of q-space where Vq-k 
is relatively small i.e., the denominator is close to one. 
For simplicity, we just set it equal to one. We then obtain 



/fe = 1 - a 



f 




1 4.^/, ''^ 


Vq — Vk Vq-k 



(19) 



This equation is much simpler to solve because the de- 
composition of Vq^k into the eigenfunctions of the square 
lattice involves only four eigenfunctions: 



Vk-q = VkVq + VkVq + l^k^q + VkVq , 



(20) 



where 



Vk = ^(cosfc^ +coBky) ; Vk = ^(cosfcj, - cosfcy) ; 

Dk = - {sinkx + smky) ; Vk = -(sin/c^ - sinky) . (21) 

We now choose a general ansatz for fk consistent with 
Eq.([9|) 

fk = A + Bvl + Gvk^k+DvkVk + Evki^k (22) 
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and solve this set of self-consistent algebraic equations 
for the coefficients. We found that the coefficients C, D 
and E are equal to zero, while A and B are the solutions 
of two coupled equations 



A = 1 - 2a 



B^2a 



(Pc 



47r2 A 



Bvl 



(23) 



Introducing A = \ — 2ax, B = lax and separating real 
and imaginary parts of x by introducing x = x\ -\- ix2, 
we obtain an equivalent set of equations for xi and X2 



Xl 



X2 



47r2 



1 - 2aa;i(l 



-I) 



vl){'2ax2) 



47r2 



1 - 2aa;i(l 



1/2) + 4a2j;2(l - jy2)2 



(24) 



In terms of xi and X2, the quasiparticle Green's function 
is given by 



G{k,u)) 



1 I - 2axi{l - ul) + i2aX2{l - vl) 

2\2 



{l-2axi{l-vl)) +^a-^xl{l-vl) 



Obviously, the spectral function and hence the DOS are 
finite when X2 ^ Q- 

A simple analysis of Eq.(^J) shows that the solution 
with X2 = exists only for \uj\ > ijJmax — 2.97tV2S 
(or a < acr — 0.448). At the critical point, we obtain 
Xl = 0.43. For smaller frequencies Eq.(^) yields a so- 
lution with finite imaginary part, just as we found with 
the toy model with momentum independent 5*. The to- 
tal bandwidth is equal to 11^ = 2uJrnax ~ 6i\/2S upto 
corrections of order 0{JS) which we neglected. For uj 
only slightly below uj„iax , we have 



X2 ^ VLU^, 



(26) 



Substituting this into Eq.(p5|), we obtain that the DOS 
behaves near ujmax as 



N{uj) 



1 



1/2 



(27) 



The above results are valid only as long as one can 
neglect the magnon dispersion. We now estimate the 
range of validity of this approximation. Recall that in 
transforming Eq.(|l3|) into Eq. (p^ , we omitted the term 



Pq 
47r2 



^ik,q - k)[G{q,LU + cOg) - G{q,u;)] 



(28) 



Far from LOmax, we do not expect this term to be rele- 
vant. Near the maximum frequency, G{u!) — G{u!max) oc 
{'-^max — ti^)^^^, and dG/duj is singular. Substituting 
the form of G from Eq.(p5|) with X2 from Eq.(^) into 
Eq.(p8|) we find that the term we omitted can be ne- 
glected when IWmax - w| > ■PS^/H/{uJma_-^ - cj)2, i.e., 
when l^niax - wl > A where A = J S{il J^/Sfl'^ . At 
frequencies closer to ujmaxj the magnon dispersion is not 
negligible, and the calculation of the spectral function 
should be done using the full self-consistency equation 
Eq.(|l3|). We will proceed with this calculation in the 
next section. 



B. Coherent part of the excitation spectrum 

In this section, we study the form of the quasiparticle 
Green's function close to the top of the valence band, i.e., 
in the region \ujmax — w| < A. 

It is again instructive to consider first a toy model 
with a momentum independent interaction. Assume 
that typical value of the magnon frequency is JS with 
J/J = 0(1). We then have instead of Eq.(|l|) 



G-^{uj)^uj-t SG{lo + JS) 



(29) 



(25) G{co) 



In the vicinity of ujmax, the solution of this equation is 

2 A / 2 ((c.-..„,„,)3+A3)i/2^ 



1 



1^^ — ^max\ 



(30) 



where Umax = 2t^/S + 0{JS) and A = JS{t/ Jy/sy/^ . 
We see that there are two typical scales introduced by 
J. For \uj — LOmax\ > A, G{ijj) in Eq.(30) differs from 
that in Eq.(p^) only by small corrections. For JS < |w — 
^max\ < A, the frequency dependence of the full solution 
is different from that in Eq. (p^) , however, G(w) remains 



approximately equal to 2/li. 



Finally, at 



< 



JS, the full Green's function begin to increase, and very 



near uj^, 



we have 



G{u) 



jVs 



t 



(31) 



We see that very near ujmax, the Green's function has 
a conventional pole with the residue Z = J^/S/i. This 
implies that around LUmax, there should exist coherent 
fermionic excitations. 

We now proceed with the solution of the actual 
self-consistency equation with a momentum-dependent 
'i'{k,q — k). Inspired by the solution of the toy model, 
we assume that there exists a frequency, tOmax for which 
G~^{k,uJmax) = at some k = k^, and which differs from 
the previously found onset frequency only by an amount 
of 0{JS). We will not be able to fully verify this assump- 
tion analytically as it would require us to find a solution 
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of Eq.(13) for all k and uj ~ Wmax which we cannot do. 
However, we will later verify this assumption in our nu- 
merical studies. We also assume and then verify that 
^0 = ('''/2, 7r/2), and that near k = kg and uj = ujmax, 
the excitations arc mostly coherent, and the quasiparticle 
Green's function has the form 



G{k,Lu) = 



Z 



^ - ^rnax + Ek - ij{u} - UJ„iax)^Q{i^ - UJrnax) 



Here Zk is the quasiparticle residue, 7 is the damping 
coefficient, Q{x) = 1(0) if a; < (a; > 0), and the hole 
excitation spectrum has the form 



{k±-ko)^ , (Mi^kof 



2m I 



2m\\ 



(33) 



where fc_L, k\\ are the momenta along the boundary of the 
MBZ and along the zone diagonal, respectively. 

In addition, as we discussed above, we introduce an up- 
per cutoff (jc < 1 hi the integration over the magnon mo- 
mentum, and restrict with an expansion of the magnon 
energy upto linear order in q. We recall that physi- 
cally, the presence of this cutoff reflects the experimen- 
tal fact that the zone-boundary magnons cease to exist 
as well-defined quasiparticles and therefore effectively do 
not contribute to the self-energy of the valence fermions. 
We will see that the quasiparticle residue Z scales as 
(qc)^^"^^, but the effective masses are independent of qc- 

We now substitute the coherent ansatz for G(fc, lo) into 
the self-consistency equation Eq. (^) . Expanding around 
fco and LOmax and using the fact that G^^{kQ,ujmax) = 0, 
we obtain self-consistent solutions for the quasiparticle 
residue, the quasiparticle spectrum and the damping co- 
efficient. Consider first the quasiparticle residue. Setting 
k = kg and expanding the r.h.s. of the self-consistency 
equation Eq. (03) to linear order in LOmax — w we obtain 



I - Z 
Z 



d?q 



4.2*(fco,9) 



Eko+,? 



(34) 



where the integration runs upto qc- Since qc <C 1, we 
can expand the two terms in the denominator to linear 
order in q. As Ug (x q and E^^^g cc q^, the first term 
is dominant. Performing the integration with only ujq in 
the denominator, we obtain 



Z 



V2PZ' 



nJ^S 



(35) 



In the limit J^/S/t <^ 1, the term linear in Z can be 
neglected and we find 



Z = 



JVS /7r\/2\i/2 



t 



^ Qc ^ 



(36) 



We see that Z scales linearly with J^/S/t as in our toy 
model. This dependence was also obtained in earlier 



studies [g. It was however noticed in Ref. 1^ that the 
linear dependence exists only for very small J /t. These 
authors argued that for moderate J/t, Z ^ (J/i)-^/-^. We 
also found deviations from the linear behavior already for 
moderately small J/i, however, we did not find a square 
root dependence for intermediate J/t. A plot of Z versus 
J/t is presented in Figs. || and 0. 

Next, we calculate the quasiparticle damping coeffi- 
cient 7. For this we again set k = ko, neglect Ekg+q com- 



(32) pared to ujq, but do not expand in lu 



nax- However, 

since we are interested in small deviations from LOmax we 
can neglect the damping term on the r.h.s. of Eq.(|T^) 
compared to w — tUmax- The r.h.s. of the self-consistency 
equation then takes the form 



^■^ik,,q)- 
47r^ UJ 



(37) 



Clearly, for lo > uJmax, the denominator is positive and 
the integral does not contain an imaginary part. For 
LU < LUmax, howcver, the integrand has a pole ai lo — 
^max — Wg. Integrating around the pole, we obtain a 
finite imaginary part which in 2D scales as {uj — oomax)'^ ■ 
After performing the explicit calculations, we obtain 



t'Z' 



(2^)2 J3 



(38) 



The same result was obtained earlier by Kane, Lee and 
Read Note in passing that in contrast to a recent 
claim in Ref. jl^], we did not find a missing factor of 2 
in their formula. Substituting the expression for Z into 
Eq. (^8|) , we finally obtain 
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1 7rV2 



4J5 qc 



(39) 



Comparing now the damping term with the term linear 
in frequency, we find that the fermionic excitations are 
weakly damped for Ek = oomax — io < AJS{qc/TT^/2). For 
small Qc, this condition is satisfied only in a small region 
around fcg. For example, for qc = \/t^/2 and ~ 4JS', 
the fermionic excitations arc only weakly damped for 
|fc — /cqI < 0.75 which constitutes only a small fraction 
of the MBZ. Away from this region, the damping term 
is dominant, and the spectral function should possess a 
broad maximum around 00 = Ek rather than a sharp 
quasiparticle peak. This is in full agreement with the 
data [|30y3l[ which show that the spectral function mea- 
sured in photoemission experiments possesses a clearly 
distinguishable quasiparticle peak only in the vicinity of 

We now proceed with the calculation of the effective 
masses. For this we set lo = Lo^ax and expand in the 
magnon momentum. We restrict ourselves to the strong 
coupling limit J^/S <C t and neglect the bare dispersion, 
which in this limit is completely overshadowed by the 
self-energy correction. We then obtain 
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Z 



4^ 



X |G(Wma:E + UJq,k + q) - G{LOmax + i^q, ko+ q)^^ 

^|*(fc,g) - '^{kQ,q)^G{uj„,ax + ujq,kQ + q) 
- J 0|vl'(fc,g)-vl/(fco,g)| 



X < G{u!r,iax +UJq,k + q)- G{uJmax +U!q,ko + q) A (40) 



Expanding the quasiparticle Green's function, we obtain 




GiuJ,nax +UJq,k + q) - G{uJ„iax + ^9, ^0 + (?) 



ko+q 



El 



k+q 



K + Eko+qY 



z 



where 

Ek+q 



Ek, 



{Ek„+q - Ek+qf 
{Uq + Ek„+qY 



k±qA 



0+9 



2 m 



2m II 



^11 



(41) 



(42) 



S 



Frequency (0 

FIG. 4. The schematic form of the DOS at half-filling. The 
DOS reaches a maximum at energies ^ JS, below the gap, 
then drops down at slightly larger energies, and then gradu- 
ally increases and saturates at energies which are ~ t away 
from the gap. 



coupling. Then only the first solution is physically rele- 
vant, and we finally obtain 



m|| = TOj_ ~ (AJS) 



(45) 



The expansion of "^(k.q) upto quadratic order in k and 
upto linear order in q yields 



^ik,q)-^'{ko,q) = 32PS 
xi^klqil 



HA] 

2V2q J 



(43) 



Inserting now these expressions into Eq.(^) and using 
the result for the quasiparticle residue, we find that one 
of the two contributions to the first term on the r.h.s. of 
Eq.(^ cancels out the Ek/Z term on the l.h.s. The 
remaining terms are all proportional to and there- 
fore the energy scale t drops completely out of the prob- 
lem. The only remaining scale is given by LOq, and the 
inverse effective masses are therefore proportional to the 
spin-wave velocity. Furthermore, we found that the in- 
tegrals over the magnon momentum in the remaining 
terms in Eq.(^o|) are confined to the upper limit of the 
q— integration, and all scale as q^. Therefore, the cutoff 
qc also drops out of the problem. As a result, we obtain 
universal, model-independent equations for the masses 



We see that in the limit qc ^ 1 and for J^/S/t <C 1 when 
the bare dispersion can be neglected, the effective masses 
are equal, i.e., the top of the valence band is isotropic. 
This result is an intrinsic property of the Hubbard (or 
t — J) model at strong coupling, independent of the form 
of the bare hopping. Note that the value of the masses 
is exactly the same as in the mean-field theory with t' = 
-0.5J. 

We also estimated the magnitude of the corrections to 
this universal result for the masses. We indeed found 
that as qc increases, the dispersion becomes anisotropic 
with m^ > m||. This trend is consistent with the results 
of other authors who integrated over the full magnetic 
Brillouin zone in Eq.(|l|) ||||,|lll- Finally, the form of the 
coherent part of the Green's function in Eq.(|3^) implies 
that the DOS behaves as {wmax — w)^ very near uJmax 
and reaches the value 

at UJmax — w ~ JS^ which is the largest scale where this 



form is applicable, 
scales as 



At even larger frequencies, the DOS 



N ' 



Z 



(47) 



1 



1 = 



1 = 



(44) 



(4J5'mi)2 JSm^ 
1 

The second equation yields m|j 

we obtain two solutions: m'^ — -^u ^ g> 
checked that only the first solution for mj^ can be contin- 
uously connected with the perturbative solution at weak 



= AJS, while for m,±^ 
US or I JS". We have 



and transforms into the fully incoherent DOS given by 
Eq.(p7|) at uomax — ui ^ K. This incoherent density of 
states then gradually increases with frequency and satu- 
rates at (jj <C iOmax- These results imply that the DOS 
reaches a maximum at oJmax — lu ^ JS, then drops 
down at slightly larger frequencies ~ A, and then grad- 
ually increases and passes through a broad maximum at 
Lu <ti UJmax- The bchavior of the DOS is presented in 
Fig. ^. This behavior is in agreement with the numeri- 
cal results which also find a strong coherent peak in the 
DOS at g^niax ~ a; ^ J on top of a smooth incoherent 



background |47 
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C. Vertex corrections 

Finally, we consider the effect of vertex corrections. 
We already found in Sec. |l| that these corrections do not 
introduce any new scale, but at the same time they do 
not possess a factor of 1/S' and therefore can only be ne- 
glected due to a numerical smallness. To estimate the 
magnitude of the vertex renormalization, we cornputed 
the lowest-order vertex correction shown in Fig, p with 
the full quasiparticle Green's functions from Eq.(32). We 
followed the same computational steps as before, i.e., ex- 
panded to linear order in the magnon momentum and 
integrated upto Qc- Performing these calculations, we 
obtain that at small external momenta the full vertex 
has the same functional form as the bare one and differs 
from it by a factor {1 + S) where 



6^ 



2ttJ^S 



(48) 



and 



I = dx dy 
Jo Jo 



xy 



{x + y) 



(log 2 -0.5) « 0.2 . (49) 



Substituting Eq.(|3q) for Z into Eq.([48|), we obtain that 
the term in brackets is equal to unity, i.e., S = I ~ 0.2. 
We see that the leading vertex correction accounts for 
only a 20% renormalization of the bare vertex. We did 
not explicitly compute higher-order vertex corrections, 
but our estimates show that they are likely to be pro- 
gressively smaller. We therefore estimate that our ana- 
lytical and numerical results for the dispersion obtained 
without vertex corrections are valid with an accuracy of 
about 20%. 



IV. NUMERICAL RESULTS 

We now proceed with the discussion of the full nu- 
merical solution of the self-consistency equation for the 
quasiparticle Green's function. 

A s in th e previous section, we begin by considering in 
Sec. IV A the frequency range \lo — LOmax\ S> A, in which 
the spectrum is completely incoherent. In Sec. IV B we 
then consider frequencies close to iOmax for which we ob- 
tain coherent excitations on the scale of 0{JS). We 



present the results for the dispersion of a single hole for 
different values of J /t and t' /t, as well as different cutoffs 
qc- For comparison with earlier studies we also present 
the results for the case when the magnons are considered 
as free particles. 

We will demonstrate that for t' — one recovers two 
equal effective masses only for very large t/J. However, 
after including a nearest-neighbor hopping t' = — 0.5J, 
we obtain two roughly equal masses for all values of t/ J . 
In this situation, the only effect of the decrease of J/t 
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(CO+A) ( in units of t) 

FIG. 5. The incoherent part of the hole excitation spec- 
trum for several values of J/t (solid line J/t — 0.4, dashed 
line J/t = 0.007 and dotted lines J/t = 0.0007). 



is the transfer of the spectral weight from the coherent 
to the incoherent part of the dispersion. Furthermore, 
for the experimentally relevant case J/t = 0.4, we find 
TO~i w {4:JSy\ For's = 1/2 this yields m-^ w (2J)~\ 
which is the same value that was obtained in the photoe- 
mission experiments on Sr2Cu02Cl2- 



A. Incoherent part of the excitation spectrum 



As we discussed in Sec. Ill A , for | w— Wmax | ^ A we can 
neglect the magnon dispersion on the r.h.s of Eq. (ITll) and 
consider an integral equation only in momentum space. 
Following the same argument we al so neglected the bare 
fermionic dispersion in Sec. Ill A . For our numerical 
studies, however, we kept the bare fermionic dispersion 
in order to illustrate how the DOS evolves with J/t. In 
Fig. ||we present for several values of J/t the DOS result- 
ing from the numerical solution of Eq.(^4|) for 5=1/2 
and t' = 0. We see that for intermediate J/t — 0.4 
(solid line), the DOS is asymmetric around w = with 
the density shifted towards negative frequencies. This 
indicates that the contribution from the bare dispersion 
which by itself yields a finite DOS only for negative ld is 
not negligible. With decreasing J/t the asymmetry be- 
comes weaker, until it basically vanishes for J/t — 0.007. 
This result is expected since in the limit J/t the 
bare dispersion becomes irrelevant, and we should re- 
cover a symmetric DOS. We also found that the total 
bandwidth only weakly depends on J/t and is roughly 
equal to = 6.6t. This value is only slightly la rger t han 
W = 6t which we obtained analytically in Sec. Ill A. 



B. Coherent part of the excitation spectrum 

In order to solve Eq. (pi]) for the full Green's function 
we use a discrete mesh in frequency and in fc— space. 
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magnon integration over 1/4 of the MBZ 




FIG. 6. Z as a function of J/t. The dotted and solid lines 
are the plots of Eqs.(p5[) and (p^), respectively. The dashed 
line is our analytical result with subleading corrections in Qc. 
The filled diamands are our numerical results. The integra- 
tion over magnon momenta is restricted to 1/4 of the MBZ. 



We assume that near the top of the valence band the 
quasiparticle Green's function has the form presented in 
Eq.(p2[) and obtain the onset frequency ujmax and the 
hole dispersion Ek from the conditions 



G ^{k = ko,UJ = UJrnax) = , 

G'~^(fc, W = UJrnax) Ek/Z 



(50) 



To obtain the quasiparticle residue, we compute G ^{k 
kf),uj) and use the relation 



Z 



G '^{k,U!jnax + Auj) - G ^{k, UJrnax) 



(51) 



where Auj is a small shift from the maximal frequency. 
The dispersion extracted from Eq.(|50|) is formally valid 
only in the vicinity of kg. At larger distances from fcoj 
Ek does not necessary coincide with the position of the 
maximum in the spectral function due to a strong quasi- 
particle damping. In our numerical procedure for solving 
the self-consistency equation, we relate the Green's func- 
tion at a given frequency lo to the Green's functions at 
larger u; + ujq, and progressively compute G at smaller 
and smaller u. Using this method, we cannot obtain the 
imaginary part of the full Green's function and therefore 
are unable to compare Ek extracted from Eq. (pO|) with 
the position of the maximum in the spectral function. 
We just assume without proof that at least not too far 
from fcoi Ek and the peak position roughly coincide. 

We first present in Figs. ^ and ^ our results for the 
quasiparticle residue at fc = fco = (7r/2,7r/2) as a func- 
tion of J/t. We have chosen two values oi Qc- a smaller 
one Qc — \/27r/16 and a larger one for which the integra- 
tion over the magnon momenta runs over 1 /4 of the MBZ. 
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FIG. 7. The same as in Fig. | for = 2^/271/16. 



The squares in these figures represent our numerical re- 
sults, the dotted line is our analytical formula, Eq. (^ , 
obtained to leading order in and the dashed line in- 
corporates subleading corrections in qc. We see that for 
smaller g^j the agreement between the numerical data 
and the results to leading order in q^ is rather good. For 
larger q^, subleading corrections are more relevant. In 
both cases, however, the quasiparticle residue is substan- 
tially reduced from its value Z — \ for free fermions 
already for moderate Jjt. We also see that the linear 
dependence exists only for very small Jjt (see linear fit 
in Fig. I). 

In Fig. H we present the results for the ratio of the 
masses as a function of i/J for i' = and t' = —0.5 J, 
respectively. In both cases, the integration over the 
magnon momenta runs over 1/4 of the MBZ. We see that 
for both values of t' , the ratio of the masses approaches 
one in the limit t/ J 00. This is in full agreement with 
our analytical results. We also see, however, that for 
t' — 0, one needs very large, unphysical values of t/J to 
recover the limiting behavior. For t' = —0.5 J, the ratio 
of the masses is equal to one already at the mean-field 
level, and our results demonstrate that the ratio remains 
roughly equal to one for all values of t/J including the 
experimentally relevant t/J = 2.5. In order to see the ef- 
fect of the J/t ratio on the whole fermionic dispersion, we 
present the results for E{k) for t' = and two different 
values of J/t in Fig. ||. We clearly see that the variation of 
J/t mainly affects the dispersion around (0,7r). The ex- 
citation energy in this region increases with t/J, i.e., the 
dip in the dispersion becomes deeper, which immediately 
leads to a decrease in the ratio of the effective masses. 
At the same time, the overall bandwidth only slightly in- 
creases with decreasing J/t. In Fig. |l^ we compare the 
results for the fermionic dispersion for J/t — 0.4 and for 
two values of t',t' = and t' = —0.4 J (here the magnon 
integration runs over 1/4 of the MBZ). We see that for 
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a) 



magnon integration over 1/4 of tiie MBZ 



b) 
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FIG. 8. The ratio of the effective masses m± and my as 
a function of t/J for (a) t' = and (6) t' = -0.5J. The 
integration over the magnon momenta runs over 1/4 of the 
MBZ. 



t' = 0, the dispersion is rather anisotropic and incon- 
sistent with the experimental data |30|j3l] |. On the con- 
trary, for t' ~ —OAJ, not only the masses are equal, but 
also the energies at (0,0) and (0,7r) are nearly equal to 
each other, and the bandwidth for coherent excitations is 
about 3 J. All three of these results are in full agreement 
with the data |]30| , ^ . We also found that the energy at 
(0,7r/2) is about half of that at (0,0) which agrees with 
the most recent data by LaRosa et al The results 

for t' = — 0.5J are similar to those for t' — —OAJ and 
are presented in Fig. ^ In this figure we compare the 
dispersion for t' = —0.5 J for two different ranges of inte- 
gration over the magnon momentum. We see that while 
the integration over 1 /4 of the MBZ yields a dispersion 
roughly consistent with the data, the integration over the 
full MBZ yields a highly anisotropic dispersion. However, 
one can increase t' even further and reduce the energy at 
(0,7r) thus making the dispersion near (7r/2,7r/2) more 
isotropic even for the integration over the full MBZ. We 
illustrate this in Fig. |lj where we present the results for 
the fermionic dispersion for t' = — J and for the integra- 
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FIG. 9. The fermionic dispersions for t' — 0.0 and two 
different values for J/t (solid line J/t = 0.4, dotted line 
J/t = 0.3). The magnon integration is restricted to 1/4 of 
the MBZ. 
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FIG. 10. The fermionic dispersions for J/t — 0.4 and 
two different values for t' (solid line t' — 0.0, dotted line 
t' — —0.4 J). The magnon integration is restricted to 1/4 of 
the MBZ. 



tion over the full MBZ. We see, however, that the overall 
bandwidth is still larger than in the experiments. We 
therefore conclude that if the integration over magnon 
momentum runs over the full MBZ (which implies that 
magnons are treated as free particles), the dispersion is 
inconsistent with the data for all reasonable values of t'. 
In this situation, to account for the data one has to adjust 
the hopping to even further neighbors. 

For completeness, we also present several results for the 
integration over the full MBZ in the conventional t — J 
model without the three-cite term. This corresponds to 
neglecting the bare fermionic dispersion in Eq.(|ll|). In 
Fig. |l^ we present the results for the excitation energy 
for <' = and J/t = 0.4. This form of the dispersion is 
in very good agreement with the results of earlier stud- 
ies As in previous studies, we found that the 
effective mass along the zone diagonal is roughly 7 times 
smaller than the mass along the boundary of the MBZ. 
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FIG. 11. The fermionic dispersions for J/t = 0.4 and 
/' = — 0.5J and two different ranges of integration. The solid 
and dotted lines are the results for the integration over the 
full MBZ and 1/4 of the MBZ, respectively. 
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FIG. 12. The fermionic dispersions for J/t = 0.4 and 
/' — ~,J. The integration over magnon momenta extends 
over the full MBZ. 



In Fig. ^ we present the results for the evolution of the 
dispersion with t' . We see that as t' increases, the ef- 
fective mass along (0, tt) direction gets smaller. How- 
ever, a rather large \t'\ is needed to reproduce two equal 
masses. Moreover, for equal masses, the overall band- 
width is about two times smaller than in the experiments. 
We see again that without restricting the integration over 
magnon momentum, one needs to add and to fine tune 
the hopping parameters to even further neighbors to re- 
produce the experimental data. 



V. SUMMARY 

We now summarize the results of our studies. We con- 
sidered in this paper the dispersion of a single hole in- 
jected into a quantum antiferromagnet. We applied a 
spin-density- wave formalism extended to large number of 
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FIG. 13. The fermionic dispersion in the t ~ J model for 
J/t = 0.4 and /' = 0. The integration over magnon momenta 
runs over the full MBZ. 



orbitals n = 25', and obtained an integral equation for the 
full quasiparticle Green's function in the self-consistent 
"non-crossing" Born approximation. At 5 = oo, the 
mean-field theory is exact. At finite S, we found that 
the self-energy correction to the mean-field formula for 
G{k,uj) scales as t/J\/S, and for large t/J, relevant to 
experiments, is small only in the unphysical limit of a 
very large spin. We found that for t/J\fS S> 1, the 
bare fermionic dispersion is completely overshadowed by 
the self-energy corrections. In this case, the quasiparti- 
cle Green's function contains a broad incoherent contin- 
uum which extends over a frequency range of ~ 6t. In 
addition, there exists a narrow region of width 0{JS) 
below the top of the valence band, where the excitations 
are mostly coherent, though with a small quasiparticle 
residue Z ^ J^/S/t. The top of the valence band is 
located at (7r/2,7r/2). 

We found that the form of the fermionic dispersion, 
and, in particular, the ratio of the effective masses 
near (7r/2,7r/2) strongly depend on the assumptions one 
makes for the form of the magnon propagator. For free 
magnons, the integration over magnon momenta in the 
self-energy runs over the whole MBZ. In this case, we 
found, in agreement with earlier studies, that the disper- 
sion around (7r/2, 7r/2) is anisotropic with a much smaller 
mass along the zone diagonal. This result holds even if 
the bare dispersion contains a sizable /' term. 

We, however, argued in the paper that the two-magnon 
Raman scattering [ p3| as well as neutron scattering ex- 
periments strongly suggest that the zone boundary 
magnons are not free particles since a substantial portion 
of their spectral weight is transformed into an incoherent 
background. Most probably, this transformation is due 
to a strong magnon-phonon interaction. In this situation, 
only magnons with small momenta, for which the interac- 
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FIG. 14. The fermionic dispersion within the t — J model 
for J/t — 0.4 and for two different values of t' (dotted line 
t' — 0.4J, solid line t' = —J). The integration over magnon 
momenta runs over the full MBZ. 



tion with phonons is necessary small, actually contribute 
to the self-energy. We modeled this effect by introducing 
a cutoff Qc in the integration over magnon momenta. We 
found analytically that for small qc, the strong coupling 
solution for the Green's function is universal, and both 
of the effective masses are equal to {4JS)~^. We further 
studied numerically the shape of the dispersion for in- 
termediate J/t and found that for t' ~ —0.5 J the ratio 
of the masses remains roughly equal to one for basically 
all values of J/t. This particular value for t' is obtained 
from a comparison of the low energy excitations in the 
underlying three-band model and the effective one-band 
Hubbard model |Q. We computed the full fermionic 
dispersion for J/t = 0.4 relevant for Sr2Cu02Cl2, and 
t' = — 0.4J and found that not only the masses are both 
equal to (2J)^^, but also the energies at (0,0) and (0,7r) 
are equal, the energy at (0, tt/2) is about half of that at 
(0,0), and the bandwidth for the coherent excitations is 
around 3 J. All of these results are in full agreement with 
the experimental data by LaRosa et al and also by 
Wells et al (except for a slightly smaller bandwidth 
and larger energy at (0,7r)) Finally, we computed the 
damping of the coherent fermionic excitations and found 
that it is small only in a narrow range around (7r/2, 7r/2). 
Away from the vicinity of (7r/2, 7r/2), the excitations are 
overdamped, and the spectral function possesses a broad 
maximum rather than a sharp quasiparticle peak. This 
last feature was also reported in the photoemission ex- 
periments. 

One of the goals of the present paper was to demon- 
strate that the experimental data for Sr2Cu02Cl2 can 
be described without introducing a spin-charge separa- 
tion. In this respect, we predict that the data should 
not change much if the experiments are performed well 
below Tc though some anisotropy of the masses is possi- 



ble because the spin damping decreases with decreasing 
T and hence qc should become larger. This prediction 
is contrary to the one derived from a model with spin- 
charge separation [Q. In this last case, it was suggested 
that the minimal model with t' — Q already accounts for 
the key experimental features, and that well below Tc, 
spinous and holons are confined such that one should re- 
cover a strong anisotropy of the masses, similar to that 
in Fig. |l| 

A final remark. Though the point of departure of our 
analysis is very different from the one in the scenario 
based on spin-charge separation ||3^], in many respects 
there exists a striking similarity between the results ob- 
tained in both approaches. First, we found that the ex- 
citations are mostly incoherent, and the bandwidth of 
incoherent excitations is several t. Second, we obtained 
that the dispersing excitations observed in photoemission 
measurements exists upto an energy scale which is given 
by J rather than by t. Both of these results are in full 
agreement with the results obtained by Laughlin in the 
framework of spin-charge separation. However, contrary 
to Laughlin, we did find a conventional Fermi-liquid pole 
in G{k,uj) near (7r/2,7r/2). The quasiparticle residue of 
the coherent excitations is small in the strong coupling 
limit and vanishes when J/t 0. In view of these results, 
we suspect that spinons and holons are actually confined 
even above the Neel temperature, but the confinement is 
weak near {tt/2, tt/2) and disappears when J/t ^ 0. A 
detailed study of this confinement is clearly called for. 
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